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Abstract 

In this paper we determine the stretch factor of the Li-Delaunay and Loo-Delaunay 
triangulations, and we show that this stretch is \/4 + 2\/2 2.61. Between any two 
points X, y of such triangulations, we construct a path whose length is no more than 



\ V 4 + 2a/2 times the Euclidean distance between x and y, and this bound is best possible. 

CsJ ■ This definitively improves the 25-year old bound of -^/lO by Chew (SoCG '86). 

To the best of our knowledge, this is the first time the stretch factor of the well-studied 
Lp-Delaunay triangulations, for any real p > 1, is determined exactly. 
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1 Introduction 



X: 
\^ 

■ Given a set of points P on the plane, the Delaunay triangulation for P is a spanning subgraph 

of the Euclidean graph on P that is the dual of the Voronoi' diagram of P. The Delaunay 
triangulation is a fundamental structure with many applications in computational geometry 
and other areas of Computer Science. In some applications (including on-line routing [B M04] ) . 
the triangulation is used as a spanner, defined as a spanning subgraph in which the distance 
between any pair of points is no more than a constant multiplicative ratio of the Euclidean 
distance between the points. The constant ratio is typically referred to as the stretch factor of 
the spanner. While Delaunay triangulations have been studied extensively, obtaining a tight 
bound on its stretch factor has been elusive even after decades of attempts. 

In the mid-1980s, it was not known whether Delaunay triangulations were spanners at all. 
In order to gain an understanding of the spanning properties of Delaunay triangulations, Chew 
considered related, "easier" structures. In his seminal 1986 paper [Che86] . he proved that an 
Li-Delaunay triangulation — the dual of the Voronoi diagram of P based on the Li-metric 
rather than the L2-metric — has a stretch factor bounded by \/T0- Chew then continued on 
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Table 1: Key stretch factor upper bounds (optimal values are bold). 



and showed that the a TD-Delaunay triangulation — the dual of a Voronoi diagram defined 
using a Triangular Distance, a distance function not based on a circle (L2-metric) or a square 
(Li-metric) but an equilateral triangle — has a stretch factor of 2 |Che89) . 

Finally, Dobkin et al. |DFS87j succeeded in showing that the (classical, L2-metric) Delau- 
nay triangulation of P is a spanner as well. The bound on the stretch factor they obtained 
was subsequently improved by Keil and Gutwin |KG92| as shown in Table [H The bound by 
Keil and Gutwin stood unchallenged for many years until very recently when Xia improved 
the bound to below 2 jXiall] . 

While progress has been made, none of the techniques developed so far lead to a tight 
bound on the spanning ratio. There has been some progress recently on the lower bound side. 
The trivial lower bound of 7r/2 ~ 1.5707 has recently been improved to 1.5846 [ BDL^ll] and 
then to 1.5932 [XZllj . 

While much effort has been made on understanding the stretch factor of Delaunay tri- 
angulations, little has been done on the Lp-Delaunay triangulations for p ^ 2. Lee and 
Wong |LW80) show that Li-,Loo-Delaunay triangulations have applications in scheduling prob- 
lems for 2-dimensional storage, and how to construct, for all real p > 1, Voronoi diagrams in 
the Lp-metric in 0(?7-logn) time |Lee80| . Delaunay triangulations based on arbitrary convex 
distance functions have been studied in [BCCSOS] , which shows that such geometric graphs 
are indeed plane graphs and spanners whose stretch factor depends only on the shape of the 
convex body. However, due to the general approach, the bounds on the stretch factor remain 
loose. For instance the bounds they obtain for L2-Delaunay triangulations are > 24. 

The general picture is that, in spite of much effort, with the exception of the triangular 
distance the exact value of the stretch factor of Delaunay triangulations based on any con- 
vex function is unknown. In particular, the stretch factor of L^-Delaunay triangulations is 
unknown for each p > 1. 



Our contributions. We show that the exact stretch factor of Li-Delaunay triangulations 
and Loo-Delaunay triangulations is \/4 + ~ 2.61, ultimately improving the 3.16 bound 
of Chew |Che86| . 

Technically, we use rectangular coordinates to prove the upper bound. We show that 
the distance in the triangulation between any source-destination pair of points lying on the 
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border of a horizontal rectangle of length x and of depth y < x is no more than (1 + \/2)x + y. 
The stretch factor bound then simply follows. In our proof, we construct the route from the 
source to the destination by maintaining two possible short paths, until we reach some special 
point (called inductive point) where we can apply our main inductive hypothesis. 

Despite the technical aspect of our contribution, we believe that our proof techniques 
may give insights into determining the stretch factor of other convex distance based Delau- 
nay triangulations. For example, let Pk denote the convex distance function defined by a 
regular /c-gon. We observe that the stretch factor of P^-Delaunay triangulations is known 
for k = 3,4 since P3 is the triangular distance function of |Che89j . and P4 is nothing else 
than the Loo-metric. Determining the stretch factor of P^-Delaunay triangulations for larger k 
would undoubtedly be an important step towards understanding the stretch factor of classical 
Delaunay triangulations. 

2 Preliminaries 

Given a set P of points in the two-dimensional Euclidean space, the Euclidean graph £ is the 
complete weighted graph embedded in the plane whose nodes are identified with the points. 
We assume a Cartesian coordinate system is associated with the Euclidean space and thus 
every point can be specified with x and y coordinates. For every pair of nodes u and w, the 
edge (n, w) represents the segment [uw] and its weight is the edge length defined in Euclidean 
distance: d2{u,w) = y/ dx{u, wY + dy{u, wY where dx{u,w) (resp. dy{u,w)) is the difference 
between the x (resp. y) coordinates of u and w. 

We say that a subgraph H oi a graph G is a t-spanner of G if for any pair of vertices 
u, V of G, the distance between u and t> in if is at most t times the distance between u and 
V in G; the constant t is referred to as the stretch factor of H (with respect to G). if is a 
t-spanner (or spanner for some t constant) if it is a t-spanner of £. 

In our paper, we deal with the construction of spanners based on Delaunay triangulations. 
As we saw in the introduction, the Li-Delaunay triangulation is the dual of the Voronoi' 
diagram based on the Li-metric di{u,w) = dx{u,w) + dy{u,w). A property of the Li- 
Delaunay triangulations, actually shared by all Lp-Delaunay triangulations, is that all their 
triangles can be defined in terms of empty circumscribed convex bodies (squares for Li or 
Loo and circles for L2). More precisely, let a square in the plane be a square whose sides 
are parallel to the x and y axis and let a tipped square be a square tipped at 45°. For every 
pair of points u,v G P, {u,v) is an edge in the Li-Delaunay triangulation of P iff there is 
a tipped square that has u and v on its boundary and contains no point of P in its interior 
(cf. |Che89) ). 

If a square with sides parallel to the x and y axes, rather than a tipped square, is used 
in this definition then a different triangulation is defined; it corresponds to the dual of the 
Voronoi diagram based on the Loo-metric doc{u,w) = max{dx{u,w),dy{u,w)}. We refer to 
this triangulation as the Loo-Delaunay triangulation. This triangulation is nothing more 
than the Li-Delaunay triangulation of the set of points P after rotating all the points by 45° 
around the origin. Therefore Chew's bound of \/T0 on the stretch factor of the Li-Delaunay 
triangulation ( |Che86] ) applies to Loo-Delaunay triangulations as well. In the remainder of 
this paper, we will be referring to Loo-Delaunay (rather than Li) triangulations because we 
will be (mostly) using the Loo-metric and squares, rather than tipped squares. 
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One issue with Delaunay triangulations is that there might not be a unique triangulation 
of a given set of points P. To insure uniqueness and keep our arguments simple, we make 
the usual assumption that the points in P are in general position, which for us means that 
no four points lie on the boundary of a square and no two points share the same abscissa or 
the same ordinate. 

We end this section by giving a lower bound on the stretch factor of Loo-Delaunay trian- 
gulations. 

Proposition 1 For every e > 0, there exists a set of points P in the plane such that the 
Loo -Delaunay triangulation on P has stretch factor at least 

\[a + 2V^-£ . 

This lower bound applies, of course, to Li-Delaunay triangulations as well. The proof of 
this proposition relies on the example of Figure [TJ 

Proof. Given 5 > 0, we define the set of points P as follows. Let point a be the origin 
and let points b, ci, and C2 have coordinates (1, — 1), {6, \/2 — 26), and {1 — 6,1 — 26), 
respectively. Additional k = ^j"^^ — i points are placed on line segment [aci] and another 
k on line segment [C26] in such a way that the difference in y coordinates between successive 
points on a segment is 6, as shown in Figures [TJ (W.l.o.g. assume that and thus k is 
an integer so that this can be done.) Let a = Po,Pi,P2,P3, ■ ■ ■ iVkiVk+i = ci be the labels, 
in order as they appear when moving from a to ci, of the points on segment [aci] and let 
C2 = 90) 91) 92) 93) • • • ) 9fc+i = hhe the labels, in order as they appear when moving from C2 to 
b, of the points on segment [C26], as illustrated in Figure [TJ 

Consider the square 5*1 of side length 1 — 6 and having a and pi on its west (left) and 
north sides, respectively (see Figure[l]o)). Since doo{a, C2) = dx{a, C2) = 1 — 6 and doo{pi, C2) = 
dy{pi,C2) = 1 — 6, point C2 is exactly the southeast vertex of square 5i. By symmetry, it 

1 - .5 




a) b) 



Figure 1: a) An Loo-Delaunay triangulation with stretch factor arbitrarily close to \/4 + 2\^. 
h) A closer look at the first few faces of this triangulation. 
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follows that for every i = 0,1,2, ... ,k, if Si is the square of side length 1 — 6 with pi and 
Pi+i on its west and north sides, then point qi is exactly the southeast vertex of Si. This 
means that all points qj with j ^ i as well as all points pj with j i,i + 1 must lie outside 
Si. Therefore, for every i = 0,1,2, k, points pi, Pi+i, and qi define a triangle in the 
-Loo-Delaunay triangulation T on P. A similar argument shows that the path qQ,qi, . . . , qk+i 
is in triangulation T as well. The triangulation T is illustrated in Figure [1^) . 

Having defined the set of points P and described its Loo-Delaunay triangulation T, we 
now analyze the stretch factor of T. A shortest path from a to 6 in T is, for example, 
a,pi,p2, ■ ■ ■ ,Pk,ci,b. The length of this path is 

d2ia,ci) + d2{ci,b) = sJdx{a,ciY + dy{a,ci)'^ + ^Jdxici,b)'^ + dy{ci,by 



'{V2 - 5)2 + ^2 + ^(1 _ 5)2 + (1 _ 25)2 

which tends to 2\/2 as 6 tends to 0. The Euclidean distance between a and b is: 



d2{a, b) = y 12 + {V2- 1)2 = ^4 - 2V2 . 

Therefore, it is possible to construct a Loo-Delaunay triangulation whose stretch factor is 
arbirtrarily close to: 

2^/2 



^ _ - V 4 + 2^/2 . 

V4-2V2 



□ 



3 Main result 

In this section we obtain a tight upper bound on the stretch factor of an Loo-Delaunay 
triangulation. It follows from this key theorem: 

Theorem 2 Let T he the Loo-Delaunay triangulation on a set of points P in the plane and 
let a and b he any two points of P. If x = doo{a,h) = m.ax.{dx{a,b),dy{a,b)} and y = 
nim{dxia,b), dy{a,b)} then 

dT{a,b) < {l + V2)x + y 
where dT{a,b) denotes the distance between a and b in triangulation T. 

Corollary 3 The stretch factor of the Li- and the Loo-Delaunay triangulation on a set of 
points P is at most 

Y^4 + 2\/2 ^ 2.6131259... 

Proof. By Theorem [21 an upper-bound of the stretch factor of an Loo-Delaunay triangulation 
is the maximum of the function 

{l + V2)x + y 
V^x2 + y2 
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over values x and y such that < y < x. The maximum is reached when x and y satisfy 

y/x = 1 + -v/2 5 and the maximum is equal to ^1 + (1 + ^/2)2 = A + 2^2. As L^- 
and Loo-Delaunay triangulations have the same stretch factor, this result also holds for 
-Li-Delaunay triangulations. □ 



To prove Theorem [2l we will construct a bounded length path in T between two arbitrary 
points a and b of P. To simplify the notation and the discussion, we assume that point a 
has coordinates (0, 0) and point h has coordinates (x, y) with < y < x. The segment [ah] 
divides the Euclidean plane into two half-planes; a point in the same half-plane as point (0, 1) 
is said to be above segment [ah] , otherwise it is helow. Let Ti , r2 , T3 , . . . , be the sequence 
of triangles of triangulation T that line segment [ah] intersects when moving from a to h. Let 
hi and li be the nodes of Ti other than a, with hi lying above segment [ab] and li and lying 
below. Every triangle Tj, for 1 < i < k, intersects line segment [ab] twice; let hi and li be the 
endpoints of the edge of Tj that intersects segment [ah] last, when moving on segment [ab] 
from a to b, with hi being above and li being below segment [ab]. Note that either hi = 
and Ti = A(/ij, /j, or li = and Tj = A(/ij_i, /ij, /j), for 1 < i < k. We also set 
ho = Iq = a, hk = h, and = h-i- For 1 < i < A:, we define Si to be the square whose sides 
pass through the three vertices of Tj (see Figure [2]); since T is an Loo-Delaunay triangulation, 
the interior of Si is devoid of points of P. We will refer to the sides of the square using 
the notation: N (north), E (east), S (south), and W (west). We will also use this notation 
to describe the position of an edge connecting two points lying on two sides a square: for 
example, a WN edge connects a point on the west and a point on the N side. We will say 
that an edge is gentle if the line segment corresponding to it in the graph embedding has a 
slope within [—1, 1]; otherwise we will say that it is steep. 

We will prove Theorem [2] by induction on the distance, using the Loo-metric, between a 
and b. Let R{a, b) be the rectangle with sides parallel to the x and y axes and with vertices 
at points a and h. If there is a point of P inside R{a,h), we will easily apply induction. 
The case when R{a, h) does not contain points of P — and in particular the points hi and li 
for < z < A: — is more difficult and we need to develop tools to handle it. The following 
Lemma describes the structure of the triangles Ti, . . . ,Tk when R{a,b) is empty. We need 
some additional terminology first though: we say that a point u is above (resp. below) R{a, b) 




Figure 2: Triangles Ti (with points a, hi, li), T2 (with points hi, /i2, and I2), and T3 (with 
points I2, /13, and I3) and associated squares 5i, ^2, and ^3. When traveling from a to 6 along 
segment [a,h], the edge that is hit when leaving Ti is [hi,li). 
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if < x„ < X and i/u > U (resp. ?/„ < 0). 

Lemma 4 // (a, 6) T and no point of P lies inside rectangle R{a,b), then point a lies on 
the W side of square Si, point b lies on the E side of square Sk, points hi, . . . ,hk all lie above 
R{a, b), and points h, . . . ,lk all lie below R{a, b). Furthermore, for any i such that 1 < i < k: 

a) Either Ti = A{hi-i,hi,li-i = li), points hi-i, hi, and li-i = li lie on the sides of Si 
in clockwise order, and hi) is a WN, WE, or NE edge in Si 

b) Or Ti = A(/ij„i = hi,li-i,li), points = hi, li, and lie on the sides of Si in 
clockwise order, and {li-i,li) is a WS, WE, or SE edge in Si. 

These properties are illustrated in Figure [2j 

Proof. Since points of P are in general position, points a, hi, and h must lie on 3 different 
sides of Si. Because segment [ab] intersects the interior of Si and since a is the origin and b 
is in cone of a, a can only lie on the W or S side of Si. If a lies on the S side then h ^ b 
would have to lie inside R{a,b), which is a contradiction. Therefore a lies on the W side of 
5i and, similarly, b lies on the E side of Sk. 

Since points hi {0 < i < k) are above segment [ab] and points li {0 < i < k) are below 
segment [ab], and because all squares Si {0 < i < k) intersect [ab], points hi,. . . ,hk all lie 
above R{a, b), and points li, . . . ,1^ all lie below R{a, b). 

The three vertices of Tj can be either hi = /ij-i, h-i, and li or /ij-i, hi, and U-i = li. 
Because points of T are in general position, the three vertices of Tj must appear on three 
different sides of Si. Finally, because and hi are above R{a,b), they cannot lie on the S 
side of Si, and because U-i and li are below R{a, b), they cannot lie on the N side of Si. 

If Ti = A(/ij_i, /ij, = li), points hi-i, hi, li must lie on the sides of Si in clockwise 
order. The only placements of points hi-i and hi on the sides of Si that satisfy all these 
constraints are as described in a). If Ti = A(/ij_i = hi,li-i,li), points hi, li, must lie 
on the sides of Si in clockwise order. Part b) lists the placements of points and li that 
satisfy the constraints. □ 

In the following definition, we define the points on which induction can be applied in the 
proof of Theorem [2j 

Definition 5 Let R{a,b) be empty. Square Sj is inductive if edge (lj,hj) is gentle. The point 
c= hj or c = Ij with the larger abcissa is the inductive point of inductive square Sj. 

The following lemma will be the key ingredient of our inductive proof of Theorem [2j 

Lemma 6 Assume that R{a, b) is empty. If no square Si, . . . ,Sk is inductive then 

dT{a,b) < {l + V2)x + y . 

Otherwise let Sj be the first inductive square in the sequence Si, S2, . . . , Sk. If hj is the 
inductive point of Sj then 

dria, hj) + ivh^ -y) < (1 + V2)xhj . 
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// Ij is the inductive point of Sj then 

dria, Ij) - yi^ < (1 + V2)xi. . 

Given an inductive point c, we can use use Lemma [6] to bound dT{cL,b) and then apply 
induction to bound c), but only if the position of point c relative to the position of point 
b is good, i.e., if x — Xc > \y — yd- If that is not the case, we will use the following Lemma: 

Lemma 7 Let R{a, b) be empty and let the coordinates of point c = hi or c = li satisfy 
< X - Xc < \y - yd- 

a) If c = hi, and thus < x — x/^. < y/^- — y, then there exists j, with i < j < k such 
that all edges in path hi, /ij+i, /ij+2j ■ ■ ■ ,hj are NE edges in their respective squares and 
X - Xh^ >yhj - y >0. 

b) If c = li, and thus < x — xi. < y — yi., then there exists j, with i < j < k such that 
all edges in path k, /j+i, /j+2, ■ ■ ■ ,lj are SE edges and x — xi^ > y — yi^ > 0. 

Proof. We only prove the case c = hj as the case c = li follows using a symmetric argument. 

We construct the path hi,hi^i,hi^2, ■ ■ ■ ,hj iteratively. If hi = /ij+i, we just continue 
building the path from /li+i. Otherwise, (/ij, /ij+i) is an edge of Tj+i which, by Lemma [U 
must be a WN, WE, or NE edge in square Sj+i. Since the S side of square Si+i is below 
R{a, b) and because x — Xh- < yh- — y, point hi cannot be on the W side of 5j+i (otherwise b 
would be inside square Sj+i). Thus {hi, /ij+i) is a NE edge. If x — x^^^-^ > yhi+i — y we stop, 
otherwise we continue the path construction from /ij+i. □ 

We can now prove the main theorem. 

Proof of Theorem [2l The proof is by induction on the distance, using the Loo-metric, 
between points of P (since P is finite there is only a finite number of distances to consider). 

Let a and b be the two points of P that are the closest points, using the Loo-metric. We 
assume w.l.o.g. that a has coordinates (0, 0) and b has coordinates (x, y) with < y < x. 
Since a and b are the closest points using the Loo-metric, the largest square having a as a 
southwest vertex and containing no points of P in its interior, which we call Sa must have b 
on its E side. Therefore (a, b) is an edge in T and dria, b) = d2{a,b) < x + y < {1 + V^)x + y. 

For the induction step, we again assume, w.l.o.g., that a has coordinates (0,0) and b has 
coordinates (x, y) with < y < x. 

Case 1: R{a,b) is not empty. 



A 



B 



C 



Figure 3: Partition of R{a, b) into three regions in Case 1 of the proof of Theorem [2l 
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We first consider tlie case when there is at least one point of P lying inside rectangle 
R{a, b). If there is a point c inside R{a, b) such that Uc < Xc and y — yd^ x — Xc (i.e., c lies in 
the region B shown in Figure [3] then we can apply induction to get dxia, c) < (1 + V^)xc + yc 
and dr(c, b) < {1 + V^){x — Xc) + y — yc and use these to obtain the desired bound for dxia, b). 

We now assume that there is no point inside region B. If there is still a point in R(a, b) 
then there must be one that is on the border of Sa, the square we defined in the basis step, or 
Sh, defined as the largest square having 6 as a northeast vertex and containing no points of P 
in its interior. W.l.o.g., we assume the former and thus there is an edge (a, c) G T such that 
either yc > Xc (i.e., c is inside region A shown in Figure [3] or y — yc > x — Xc (i.e., c is inside 
region C). Either way, dxia, c) = ^2(0, c) < Xc + yc- If c is in region A, since x — Xc > y — yc, 
by induction we also have that dxic, 6) < (1 + V^)(x — Xc) + {y — yc)- Then 

dxia^b) < dria, c) + dxic^b) 

< Xc + yc + {l + V2){x- Xc) + {y- yc) < (1 + V2)x + y 

In the second case, since x — Xc < y — yc, by induction we have that dric, b) < {1 + V2){y — 
yc) + {x — Xc)- Then, because y < x, 

dT{a,b) < dT{a,c) + dT{c,b) 

< Xc + yc + {l + V2){y - yc) + {x - Xc) < (1 + V2)x + y 

Case 2: The interior of R{a, b) is empty. 

If no square Si, S2, - - - , Sk is inductive, dxia, b) < {1 + \f2)x + y by Lemma[6l Otherwise, 
let Si be the first inductive square in the sequence and suppose that hi is the inductive point 
of Si. By Lemma[71 there is a j , i < j < k, such that hi, /ij+i, /ii+2; - - - ,hj is a path in T of 
length at most {xhj — XhJ + (y^. — yh^) and such that x — Xhj > yhj — y > 0. Since hj is closer 
to b, using the Loo-metric, than a is, we can apply induction to bound dT{hj,b). Putting all 
this together with Lemma [6l we get: 

dT{a,b) < dT{a,hi) + dT{hi,hj) + dT{hj,b) 

< (1 + V2)xh^ - {yh^ - y) + {xhj -Xh,) + {yh, - yhj) + {l + V2){x - Xhj) + {yhj - y) 

< (1 + V2)x . 

If li is the inductive point of Si, by Lemma [7| there is a j, i < j < k, such that 
li, li^i,li+2, - - - ,lj is a path in T of length at most (x/,^ - Xhi) + {yhj - yhi) and such that 
X — Xhj > y — yhj > 0. Because the position of j with respect to b is good and since Ij is closer 
to b, using the Loo-metric, than a is, we can apply induction to bound dT{lj,b). Putting all 
this together with Lemma [H we get: 

dT{a,b) < dT{a,li) + dT{li,lj) + dT{lj,b) 

< (1 + V2)xi^ + yi^ + [xi^ - xzj + {yi^ - yi^) + (1 + V2){x - xi^) + (y - yi.) 

< (1 + V2)x + y . 

□ 

What remains to be done is to prove Lemma [H To do this, we need to develop some 
further terminology and tools. Let Xi, for 1 < i < k, he the horizontal distance between point 
a and the E side of Si, respectively. We also set xq = 0. 
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Definition 8 A square Si has potential if 

dT{a,hi) + dT{a,li) + dsi{hi,li) < 4rcj 

where dsi{hi,li) is the Euclidean distance when moving from hi to li along the sides of Si, 
clockwise. 

Lemma 9 // R{a, b) is empty then Si has potential. Furthemore, for any 1 < i < k, if Si 
has potential but is not inductive then S'^+i has potential. 

Proof. If R{a,b) is empty then, by Lemma [U a lies on the W side of Si and xi is the side 
length of square Si. Also, hi lies on the N or E side of Si, and li lies on the S or E side of 
Si. Then dT{a,hi) + dT{a,li) + ds^{hi,li) is bounded by the perimeter of square 5i which 
is 4xi. 

Now assume that Si, for 1 < i < k, has potential but is not inductive. Squares Si and 
Si+i both contain points k and hi. Because Si is not inductive, edge {li,hi) must be steep 
and thus dx{li, hi) < dy{li, hi). To simplify the arguments, we assume that li is to the W of 
hi, i.e., X;- < x/i- . The case xi. > Xh. can be shown using equivalent arguments. 

Since Tj = A(/ij_i, /ij, = li) or Tj = A(/ij_i = hi,li-i,li), there has to be a side of 
Si between the sides on which li and hi lie, when moving clockwise from li to hi. Using the 
constraints on the position of hi and li within Si from Lemma H] and using assumptions that 
{li, hi) is steep and that xj- < x/i-, we deduce that U must be on the S side and hi must be 
on the N or E side of Si. 

If hi is on the N side of Si then, because xi- < x/^., hi must also be on the N side of ^j+i 
and either li is on the S side of 5j+i and 

dsr+Ahiik) - ds^{hi,li) = 2{xi+i-Xi) (1) 

as shown in Figured^) or li is on the W side of S'j+i, in which case 

ds,+iihi,li) - ds^ihi,li) < 4{xi+i-Xi) (2) 

as shown in Figure (Da). 

If hi is on the E side of Si then, because Xj+i > Xi (since {k, hi) is steep), hi must be on 
the N side of Sj+i and either li is on the S side of and inequality ([T|) holds or li is on the 
W side of Si+i and inequality ([2]) holds, as shown in Figured]:;). 

Since Si has potential, in all cases we obtain: 

dT{a,hi) + dT{a,li) + ds^+-^{hi,li) < 4x^+1. (3) 

Assume Tj+i = A{hi, /li+i, k = h+i); in other words, {hi, /ij+i) is an edge of T with hi+i 
lying somewhere on the boundary of Si^i beween hi and li, when moving clockwise from hi 
to li. By the triangular inequality, d2{hi, /ij+i) < dsi^^{hi, /ij+i) and we have that: 

dT{a,hi+i) + dT{a,li+i) + t^s.+i (/ii+i, ^i+i) < dT{a,hi) + dT{a,li) + d2(/ii, /ij+i) + ds^+^{hi+i,li) 

< dT{a,hi) + dT{a,li) + ds^^^{hi,li) 

< 4xi+i . 
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Figure 4: The first, second and fourth case in the proof of Lemma O In each case, the 
difference dsi+i (^i; ^i) — dsi{hi,li) is shown to be at most 4Aa;, where A^; = Xj+i — Xj. 

Thus Sj+i has potential. The argument for the case when Tj+i = A(/ij = /ij+i, /j, /j+i) is 
symmetric. □ 



Definition 10 A vertex c (hi or li) of Ti is promising in Si if it lies on the E side of Si. 

Lemma 11 If square Si has potential and c = hi or c= k is a promising point in Si then 

dxia, c) < 2xc ■ 

Proof. W.l.o.g., assume c = hi. Since hi is promising, Xc = Xh^ = Xi. Because Si has 
potential, either dT{a,hi) < 2xh. or dxiaji) + dsi{li,hi) < 2xh^. In the second case, we can 
use edge {U, hi) and the triangular inequality to obtain dxia, hi) < dxia, k) + \lihi\ < 2x/i. . □ 



Here we define the maximal high and minimal low path. 



Definition 12 

• If hj is promising in Sj, the maximal high path ending at hj is simply hj; otherwise, 
it is the path hi, /ij+i, . . . ,hj .such that /ij+i, . . . ,hj are not promising and either i = 
or hi is promising in Si. 

• If Ij is promising in Sj , the maximal low path ending at Ij is simply Ij ; otherwise, it is 
the path li, k+i, . . . ,lj such that k+i, . . . ,lj are not promising and either i = or li is 
promising in Si. 

Note that by Lemma [H all edges on the path hi, /ij+i, . . . ,hj are WN edges and thus the 
path length is bounded by {xhj - XhJ + {yt^ -yhi)- Similarly, all edges in path k, k+i, . . . ,lj 
are WS edges and the length of the path is at most {xi. — x; - ) + (y; - — yi. ) . 

We now have the tools to prove Lemma [6j 

Proof of Lemma[6l If R{a, b) is empty then, by LemmaSl b is promising. Thus, by Lemma[9] 
and Lemma mi if no square Si, . . . ,Sk is inductive then dria, 6) < 2x < (1 + \/2)x + y. 

Assume now that there is at least one inductive square in the sequence of squares 
5i, . . . , 5fc. Let Sj be the first inductive square and assume, for now, that hj is the inductive 
point in Sj. By Lemma El every square Si, for i < j, is a potential square. 
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Since {Ij, hj) is gentle, it follows that d2{lj, hj) < ^/2{xhj — xi-). Let h+i, . . . , Ij-i = Ij 
be the maximal low path ending at Ij. Note that dT{li,lj) < {xi. — a;/ J + (y;. — yi-). Either 
li = Iq = a oi li is a promising point in potential square Si] either way, by Lemma [TTl we 
have that dT{a,li) < 2xi-. Putting all this together, we get 

dT{a,hj) + {yhj - y) < dT{a,li) + dxili^lj) + d2{lj,hj) + yn^ 

< 2xi^ + (x/^, -xi.) + {yi^ - yi^) + V2{xh^ - xi.) + y/j. 

< V2xh^ + xi^ + yh^ - yi^ 

< {1 + V2)xh^ 

where the last inequality follows xi^ + yh^ — yi^ < Xhj, i.e., from the assumption that edge 
{Ij, hj) is gentle. 

If, instead, c = Ij is the inductive point in inductive square Sj, let hi, hi+i, . . . , hj^i = hj 
be the maximal high path ending at hj. Then drihi, hj) < {xh. — x^.) + {y^. — yhi)- Just as 
in the first case, we have that dj-^a, hi) < 2xfi. and 

dT{a,lj)-y^ < dT{a,hi) + dT{hi,hj) + d2{hj,lj) - yi. 

< 2xh, + {xh, -xh,) + {yhj -yh,) + V2{xi^ - xhj) - yij 

< \Pix^ + Xh^ + yh^ - yi. 

< (1 + V2)xi^ . 

where the last inequality follows from Xh + yh — Vi ^ xi., i.e., from the assumption that 
(hjjlj) is gentle. □ 



4 Conclusion and perspectives 

The -Li-Delaunay triangulation is the first type of Delaunay triangulation to be shown to 
be a spanner [Che86j . Progress on the spanning properties of the TD-Delaunay and the 
classical L2-Delaunay triangulation soon followed. In this paper, we determine the precise 
stretch factor of an Li- and Loo-Delaunay triangulation and close the problem for good. The 
techniques we develop in this paper have potential to be successfully applied to Delaunay 
triangulations defined by other regular polygons, and possibly even to the classical Delaunay 
triangulation. 

From a routing perspective, it is of interest to construct routes in geometric graphs that can 
be determined locally from a neighbor's coordinates only [BCD09] . Unfortunately, the route 
that is implicitely constructed in our proof is built using non-local decisions. It would be inter- 
esting to know whether in the Li-Delaunay triangulation a route with stretch \/4 -|- 2\/2 can 
be constructed using a local routing algorithm. For TD-Delaunay triangulations, |BFvRV12] 
showed that there is no local routing algorithm that achieves a stretch that is less than 
5/\/3 ~ 2.88, whereas the stretch factor is actually 2. We leave open the questions regard- 
ing the gap between the stretch factor of Li -Delaunay triangulations and the stretch that is 
possible using local routing. 
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